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Neutrinos propagating in media (matter and electromagnetic fields) undergo flavor and helicity oscillations, 
where helicity transitions are instigated both by electromagnetic fields and matter currents. In addition, it 
has been shown that correlations between neutrinos and antineutrinos of opposite momentum can build up 
in anisotropic media. We rederive the neutrino equations of motion in the mean-field approximation for homo¬ 
geneous yet anisotropic media, confirming previous results except for a small correction in the Majorana case. 
Furthermore, we derive the mean-field Hamiltonian induced by neutrino electromagnetic interactions. We also 
provide a phenomenological discussion of pair correlations in comparison with helicity correlations. 
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I. INTRODUCTION 

Neutrino flavor conversion in vacuum [1, 2], in matter 
[3-8], or self-induced flavor conversion in a gas of interact¬ 
ing neutrinos [9-27] provide a rich phenomenology of very 
practical experimental and astrophysical importance. The data 
leave no room for doubt that neutrinos have small but nonvan¬ 
ishing masses. One consequence is that neutrinos have small 
electromagnetic dipole and transition moments [28]. These 
lead to spin and spin-flavor oscillations in strong electromag¬ 
netic fields [29-32]. Actually, polarized matter or matter cur¬ 
rents alone instigate spin and spin-flavor transitions of mas¬ 
sive neutrinos, having effects similar to electromagnetic fields 
[33-37], 

For neutrinos streaming from a supernova core, the back¬ 
ground medium may contain currents. Moreover, the neu¬ 
trino stream itself provides an unavoidable nonisotropic back¬ 
ground. In addition, self-induced flavor conversion in an in¬ 
teracting neutrino gas requires unstable modes in flavor space 
(run-away solutions). If such solutions exist, even small per¬ 
turbations or otherwise small effects can grow exponentially. 
In this sense, it is never obvious if a seemingly small effect 
can get amplified by an instability to play an important role 
after all. Therefore, it is interesting to study if an interact¬ 
ing neutrino gas can amplify helicity conversion effects [38] 
which otherwise are very small. 

Flavor oscillations lead to correlations building up between 
neutrinos of different flavor. If a*. is the creation operator 
of a neutrino in flavor state a with a certain momentum p, 
the initially prepared system can be described by the occu¬ 
pation number («j v a 0: ). One way of looking at flavor oscil¬ 
lations is that “flavor off-diagonal” occupation numbers of 
the type (a' a ap) develop and oscillate [39—41 ]. One uni¬ 
fies these expressions in a density matrix p with components 
Pap = {a\cL a ). It evolves according to the commutator equa¬ 
tion ip = [H,p], where H is the Hamiltonian matrix, con¬ 
sisting of oscillation frequencies. For vacuum oscillations we 
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have H = M 2 /2 E, where M 2 is a matrix of squared neutrino 
masses. Similar descriptions pertain to spin and spin-flavor 
oscillation, where the indices now indicate various states of 
spin and/or flavor. 

It was recently stressed that yet another form of correla¬ 
tions, hitherto neglected in the context of neutrino propaga¬ 
tion, can build up in nonisotropic media [42^14]. If a* is 
the creation operator of a massless neutrino in mode p and 
b]_ the one for an antineutrino with opposite momentum, cor¬ 
relators of the form k p = (h_ p a p ) and k j, = (a*6L p ) will 
build up, the latter corresponding to the creation of a particle- 
antiparticle pair with vanishing total momentum. Because 
massless neutrinos and antineutrinos have opposite helicity, 
this pair has total spin 1 so that its creation requires a medium 
current transverse to p to satisfy angular-momentum conser¬ 
vation. This requirement is analogous to the case of helicity 
transitions where we also need a transverse current or mag¬ 
netic field for the same reason. Including flavor and spin de¬ 
grees of freedom expands the “pair correlations” n and K* to 
become matrices similar to p. 

To develop more intuition about the meaning of the pair cor¬ 
relations, we consider a single mode p of neutrinos and —p of 
antineutrinos. We define p p = (a p o p ) and for antineutrinos 
p p = (b_ p h- p ) involving the opposite momentum. Follow¬ 
ing the earlier literature [42-44], we unify these expressions 
in an extended density matrix 


R = 




( a p a p) ( b -p a p) \ 

<4 6 -p) ( b -p b -p)J 


(1) 


which obeys an equation of motion of the form [42-44] 

iR = [H, R]. (2) 


If the background is a current moving in the transverse direc¬ 
tion with velocity /?, the Hamiltonian matrix is found to be 


H = E 




( 3 ) 


where E = |p|. For or v T neutrinos the usual matter po¬ 
tential is V = Gpn n /\/2, where n n is the neutron density. 

This commutator equation has the same structure that one 
encounters for the evolution of any two-level system and in 
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particular for two-flavor or helicity oscillations, of course with 
a different matrix H for each case. However, what specifically 
are the two states that are being mixed by the matter current 
in the pair-correlation case? 

The answer becomes evident if one considers the evolu¬ 
tion of states rather than correlators. Our simple system is 
described by the four basis states |00), 110), 101) and 111), 
where the first entry refers to z^(p) and the second to p). 
A homogeneous background medium cannot mix states of 
different total momentum, so the single-particle states must 
evolve independently as idt 110) = (£ , +H)|10) and*<9 t |01) = 
(E — V) |01), i.e., they simply suffer the usual energy shift by 
the weak potential of the medium. This leaves us with 100) 
and 111) which both have zero momentum and therefore can 
be mixed by a homogeneous medium. The former has spin 0, 
the latter spin 1, so for the medium to mix them, it must pro¬ 
vide a transverse vector in the form of a current or a spin po¬ 
larization. If Aqq and An are the amplitudes of 100) and 111), 
respectively, we will show later that Eq. (2) corresponds to 

*(£)-(£ £)(£)• 

Therefore, it is the empty and the completely filled states that 
are being mixed and that oscillate. The true ground state of 
our system is not 100), but a suitable combination of |00) and 
111) which follows from diagonalizing the matrix in Eq. (4). 

As we have noted, any two-level system is equivalent to an 
abstract spin--) system. In two-flavor oscillations, the “spin” 
represents the two flavor states. In the pair-correlation case, 
“spin up” means “empty” and “spin down” means “full with a 
pair.” This interpretation is analogous to Anderson’s “pseudo 
spin” devised to describe Cooper pairs in the context of su¬ 
perconductivity [45]. A coherent superposition of these two 
spin states, represented in our case by the pair correlations, 
corresponds to a coherent superposition of 100) and j 11). 

In analogy to the example of superconductivity, another 
way to think about these phenomena is in terms of Bo- 
golyubov transformations of the creation and annihilation op¬ 
erators. If we think of a single momentum mode p of mixed 
neutrinos in vacuum, the operators <z„ e and a v in the flavor 
basis are rotated by a unitary transformation with mixing an¬ 
gle $ to form new operators c$a Vt . +s^a UiL and 
and similarly for the creation operators, to form a new set of 
canonically anticommuting operators, now describing neutri¬ 
nos in the mass basis. Describing flavor oscillations in terms 
of time-dependent Bogolyubov transformations can be espe¬ 
cially illuminating to understand quantum statistics in mixing 
phenomena for both bosons and fermions [46]. Pair correla¬ 
tions correspond to the same idea where the mixing is between 
a p and b_ with a mixing angle corresponding to the unitary 
transformation that diagonalizes the matrix in Eq. (4). The 
state 100) defined in the Bogolyubov-transformed basis is the 
ground state of the system and no longer oscillates into the 
new 111) state. 

The goal of our paper is two-pronged. On the one hand 
we reconsider the mean-field equations of motion for massive 
neutrinos propagating in a background medium that can con¬ 
sist of matter and neutrinos, and that is homogeneous but not 


isotropic. Besides the usual flavor oscillations in matter, the 
resulting phenomena include spin and spin-flavor oscillations 
as well as pair correlations. 

As a second goal, we provide a phenomenological discus¬ 
sion of the interpretation of the pair correlations in the context 
of neutrino oscillation problems in dense media. Ultimately, 
our community needs to develop an understanding if, from a 
practical perspective, we need to worry about pair correlations 
and helicity oscillations in the supernova context. 

The supernova environment is characterized by small neu¬ 
trino energies of at most some 200 MeV (for degenerate v e ), 
i.e., small compared to W and Z masses so that it suffices to 
describe neutrino interactions in terms of an effective current- 
current Hamiltonian. In the early Universe, where the chem¬ 
ical potentials of background particles are small, one has to 
worry about corrections from the electroweak gauge-boson 
propagators even at low temperatures [47]. The supernova 
environment, in contrast, has large densities of background 
particles and this concern is moot. 

On the mean-field level, the current of background parti¬ 
cles is a classical quantity. For example, the neutral-current 
interaction of a neutrino with neutrons is given by the Hamil¬ 
tonian density % = V2Gy where Gy is the 

Fermi constant, v is the neutrino Dirac field, P], is the left- 
handed projector, and is the neutron current. If the cur¬ 
rent is homogeneous, H = J d 3 xH is effectively a “for¬ 
ward” Hamiltonian: it couples, e.g., a' and a of equal mo¬ 
menta. Following the previous literature [40, 41, 44], the evo¬ 
lution of, e.g., the annihilator for a neutrino of mass eigen¬ 
state i in mode p is given by the Heisenberg equation of 
motion p) = [a,i(t, p), H], It is then straightfor¬ 

ward to find the equations of motion of bilinears of the form 
a] (t, p )dj (f, p), of their expectation value (a| (t, pja, (t, p)), 
of the entire matrix p, and then of the extended matrix R which 
also includes pair correlations. 

It is largely a cumbersome bookkeeping exercise to obtain, 
for neutrinos with mass, all the components of the Hamil¬ 
tonian matrix H appearing in Eq. (2) when R involves all 
components of spin and flavor. We perform this task sepa¬ 
rately for Dirac neutrinos in Sec. II, for Majorana neutrinos 
in Sec. Ill, and for Weyl neutrinos (massless two-component 
case) in Sec. IV. These derivations closely parallel the recent 
paper by Serreau and Volpe [44] and we will largely follow 
their notation to avoid confusion. In the Majorana case, we 
find a small correction, but otherwise our results agree. 

The density matrix formalism allows one to treat helicity 
oscillations induced by magnetic fields and by matter cur¬ 
rents on equal footing for both Dirac and Majorana fermions. 
We derive the mean-field Hamiltonian induced by electro¬ 
magnetic fields in Sec. V. Concerning helicity oscillations in¬ 
duced by matter currents, which we analyze in Sec. VI, our 
results coincide with those of Volpe and Serreau, and paral¬ 
lel those of Vlasenko, Fuller, and Cirigliano [36-38] as far 
as the mean-field limit is concerned. These authors have de¬ 
rived the neutrino kinetic equations starting from first field- 
theoretic principles and have carried the results beyond the 
mean-field limit to include (nonforward) collision terms, gen¬ 
eralizing previous derivations [39—41]. We note in passing 
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that one of their findings — helicity oscillations in a non¬ 
isotropic matter background — had been anticipated in sev¬ 
eral papers by Studenikin and collaborators who have worked 
out the one-to-one correspondence to the effect of electro¬ 
magnetic fields [33, 35], Of course, Vlasenko, Fuller, and 
Cirigliano also included neutrino-neutrino interactions as an 
agent of helicity conversion and carried their results beyond 
the mean-field limit. 

Pair correlations have been studied in detail in condensed 
matter and nuclear physics, as well as in the context of pair 
creation in quantum field theory. On the other hand, in neu¬ 
trino physics these concepts are less familiar. They have 
been addressed only in a handful of papers in the context 
of leptogenesis, where pair correlations have been studied 
from first principles in a series of papers by Fidler, Herra- 
nen, Kainulainen and Rahkila [48-52]. In the context of neu¬ 
trino propagation in supernovae, the only discussions so far 
appear in a series of papers by Volpe and collaborators [42- 
44]. We address phenomenological aspects of pair correla¬ 
tions in Sec. VII and compare them to helicity correlations. 

Finally, in Sec. VIII we summarize the results and present 
our conclusions. 


II. DIRAC NEUTRINO 


and which involve opposite-momentum modes. 

The kinetic equations for Eqs. (5) and (6) are obtained with 
the Heisenberg equation of motion. As we will show below, 
in the mean-field approximation, and assuming spatial homo¬ 
geneity, the Hamiltonian of charged- and neutral-current neu¬ 
trino interactions can be written in the compact form 

H mi = j d 3 xDi(t,x)Tiji/j(t,x), (7) 

where summation over repeated indices is implied. The kernel 
takes account of the background medium and is 

r a = 7 vPlV* , (8) 

where P L = (1 — 7s)/2 is the usual left-chiral projector. The 
current of background matter will be defined in Eq. (24). 

The momentum-mode decomposition of a Dirac neutrino 
field reads 

i/i(f,x)= f e lp x i/ i>s (f, p), (9) 

J p,s 

where f s denotes the phase-space integration f G? 3 p/(27r) 3 
and the summation over helicities. In the mass basis, the indi¬ 
vidual momentum modes are 


Our first goal is to derive the components of the Hamilto¬ 
nian matrix H which governs the evolution equation (2) for the 
extended density matrix R including flavor, helicity, and pair 
correlations. In this rather technical section, we begin with 
the conceptually simplest case of three neutrino flavors which 
are assumed to have Dirac masses. Therefore, helicity cor¬ 
relations involve the sterile components of the neutrino field, 
which otherwise are completely decoupled. 

A. Two-point correlators and kinetic equations 


Vi >s (t, P) = a i)S (t, p)u i>s (p) + b[ s (t, -p)v i>s (-p). (10) 

The chiral spinors u and v are given in Appendix A, and the 
creation and annihilation operators satisfy the usual equal¬ 
time anticommutation relation, 

p),a] ft (f, k )} = (27r) 3 5(p-k) l 5 ii 5 sh . (11) 

Similar relations hold for the antiparticle operators b and fot. 

As a next step, we contract the kernels Ty with the spinors 
appearing in the mean-field Hamiltonian (7), leading to the 
matrices [44] 


In the simplest approximation, one can describe the state of 
a neutrino gas in terms of one-particle distribution functions. 
They are extended to include flavor and helicity coherence ef¬ 
fects by promoting the one-particle distribution functions to 
density matrices [36, 39-41, 44]. In terms of the usual cre¬ 
ation and annihilation operators, their components are 

(2-7r) 3 <5(p—k)py iS /j(f, p) = (at h (t, +k)a i)S (f, +p)}, (5a) 

(2n) 3 S(p-\t)p ij}Sh (t, p) = (b\ s {t,-p)b jth {t, —k)), (5b) 

where i and j are flavor indices in the mass basis, and s and 
h € {+, —} denote helicities. In this convention, the density 
matrix for antineutrinos Pij tS h(t, p) for momentum p actually 
corresponds to the occupation numbers of antineutrinos with 
physical momentum p. 

This convention is necessary to combine p and p with the 
pair correlations which are defined as [43, 44] 

(27r) 3 i5(p-k)Ky >?l (f, p) = {b jt h(t, -k)a ijS (f, +p)), (6a) 

(27t) 3 <5( P —k)«t sh (t, p) = (at h (t, +p)b[ s (t, -k)), (6b) 


'ij,8h( P) = u *,»(P) r *f u j,k(p) > 

(12a) 

ij,sh(P) = ^i,s{p)^ij v j,h{~ P) i 

(12b) 

'ij,sh( p) = «i,*(-p)r i jUj ! / l (p), 

(12c) 

p) = «i,*(-p)rij^,h(-p) , 

(12d) 


in component form. We can now bring Eq. (7) to the desired 
form bilinear in the creation and annihilation operators 


ffmf 


a i, s (p) r D>h(p) a f,h(p) 

+4 s (p) r ^(p) 6 L(~p) 

+ bi ! s(-p)T^ sh {p) a j,h(p) 


+ bi,s(-P) T ij,sh(P) b j,h(~P) 


( 13 ) 


where we have omitted the time arguments to shorten the no¬ 
tation. Summation over repeated indices is implied. 

Using the Heisenberg equation of motion with this Hamil¬ 
tonian one finds the extended equation of motion iR = [H, R], 
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see also Eq. (2). The extended density matrix, Eq. (1), and the 
Hamiltonian, Eq. (3), generalize to [42] 

«-(}.:,) - h- , <»> 

where the submatrices Y\ vu = , W vv = Y vv etc. and p, k, 

etc. are 6x6 matrices in helicity and flavor space. The product 
between such matrices in the commutator is defined in the 
obvious way {A ■ B)i^ s h = Ai ntSr B n j tr h with a summation 
over repeated indices. In the following we write the matrix 
structure in the form of 2 x 2 matrices in helicity space. 


— 

ij 

-+ 



+- 

ij 

++ 


where each entry is itself a 3x3 matrix in flavor space. 

B. Hamiltonian in the mean-field approximation 

After having established the overall structure of the kinetic 
equations we now turn to the interactions contributing to neu¬ 
trino refraction in the supernova environment. In this subsec¬ 
tion we only consider charged- and neutral-current neutrino 
interactions, whereas the analysis of the electromagnetic in¬ 
teractions is postponed to Sec. V. 

1. Charged-current interaction 

We begin with charged-current (cc) interactions with back¬ 
ground charged leptons. In the low-energy limit and after a 
Fierz transformation, the usual current-current Hamiltonian 
density is 

H cc = V2 G f J2 [vcSfPLVp] $87m(1 ^ 7 5 )4] , (16) 

a,/3 

where a, /? G {e, p, r} are flavor indices. 

To obtain the neutrino mean-field Hamiltonian we replace 
the second bracket by its expectation value. In the supernova 
environment, the temperature is too low to support a substan¬ 
tial density of muons or tauons, and we use only the electron 
background. Then we find in the mass basis 

?Cf = V^Gf X; [Vi'fPLVj] [Utl^Uej] , (17) 

where U is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) 
matrix. We have introduced a linear combination of vector 
and axial-vector charged electron currents, 

!cc = evince) - CA(e7 M 7 5 e), (18) 

where cy = ca = 1. Because electrons are the only back¬ 
ground particles contributing to charged-current interactions 
and to simplify the notation, an “e” index is implied in Ig.. If 
the electrons are not polarized, the axial current vanishes and 
Ig c = ./// , the “convective” electron current. 


2 . Neutral-current interaction with matter 

The neutral-current (nc) interactions with matter are de¬ 
scribed in the mass basis by the Hamiltonian density 

n nc = V2G f Y [vilri P Lvi\ ['0/7 M ( c v r— c a7 5 )'0/] , (19) 

where / denotes electrons, protons, and neutrons. The result¬ 
ing contribution to the mean-field Hamiltonian is 

?Cf = V2G f Y [vil»PLVi][l£ c + 1$ + /£], (20) 

i 

where Ig c denotes the electron neutral current (index e im¬ 
plied), whereas the other contributions refer to protons and 
neutrons as explicitly indicated. 

These currents are defined in analogy to Eq. (18) with the 
appropriate coupling constants. For electrons, they are given 
by cy = —2 + 2 sin 2 0w (Weinberg angle 8\y) and ca = — \- 
For protons, cy = \ — 2 sin 2 0\,y, i.e., the same as for elec¬ 
trons with opposite sign, and for neutrons cy = — 2. For the 
nucleon axial vector one often uses ca = ±1.26/2 in analogy 
to their charged current. However, the strange-quark contribu¬ 
tion to the nucleon spin as well as modifications in a dense nu¬ 
clear medium leave the exact values somewhat open [53, 54]. 

In an unpolarized and electrically neutral environment, the 
axial currents disappear and the electron and proton contri¬ 
butions to the convective neutral current cancel such that in 
Eq. (20) we have Ig c +Ig+I% = where J£ is the neu¬ 

tron convective current. Neutrino refraction in such a medium 
depends only on the charged electron current and the neutral 
neutron current. 

3. Neutrino-neutrino interaction 

The most complicated interaction is the neutral-current 
neutrino-neutrino one. It is described in the mass basis by 
the Hamiltonian density 

U vv = 77| G f [^OV-Pz,^] . (21) 

V ij 

To obtain the mean-field Hamiltonian bilinear in the neutrino 
fields we need to replace products of two of the four neutrino 
fields in this expression by their expectation value. 

The only combinations that do not violate lepton number 
are of the type (Pji/j) and (viVj), where i and j can be equal 
or different. We denote the corresponding mean field as 

= {vj'fPLVi) ■ ( 22 ) 

To simplify notation we avoid an explicit “neutrino” and “nc” 
index, i.e., expressions of the type always refer to the neu¬ 
tral neutrino current for the mass states i and j. An explicit 
expression in terms of the density matrices and pair correla¬ 
tors will be given in Eq. (27) below. 

For the i = j contractions it is sufficient to take the expecta¬ 
tion value of one of the square brackets in Eq. (21), leading to 
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the mean-field Hamiltonian \/2Gf [Pi7ii P L^i\ Ijj■ For 
the i ^ j contractions we use the Fierz identity to rewrite the 
Hamiltonian as [£^ 7 j\[vjAP lA in Eq. (21), leading to 
the contribution \/2 Gf Aj- Altogether, we 

find 


where p is a unit vector in the momentum direction and the 
complex polarization vector e spans the plane orthogonal to p 
(see Appendix A for more details). We also use cf> to denote 
the polar angle of p in spherical coordinates. To lowest order 
in rrij, the spinor contractions are then found to be 


nz = v^Gf £ [d^Pluj] [7g + Sij £ fc J&] (23) 

ij 

for the neutrino-neutrino mean-field Hamiltonian. 


C. Components of the Hamiltonian matrix H 


Adding up Eqs. (17), (20), and (23) we find the overall 
mean-field current 

V?i = V2 Gf [u} e ic c u e j + Sij(+1$ + K) 


+ l ij + &ij Efe J fefe] ■ (24) 

The spinor contractions defined in Eq. (12) lead to the com¬ 
ponents of the Hamiltonian matrix H of the form 

H ZsH = (7 + SshSijEi, (25a) 

H tsh = (iMtshVt, , (25b) 

H Zsh = , (25c) 

H Zh = (7 M\l ah Vtj ~ SshS^E, , (25d) 


where Ej = (p 2 + rnj) 2 is the neutrino energy, and we 
have identified W vv = r vu , H" ! ' = Y vv , etc. We have used 


the compact notation 

( 7 ^ p L)ij, sh = u-iA+P)l» p LUjA+P ), (26a) 

( 7 » p L)ij,sh = UiA+P)ln p LVj,h(-P) > (26b) 

( 7 H P L)ij,sh = Vi,s{-Ph» P LUj,h{+P) : (26c) 

(7a i P L)ij,sh = Vi lS {-P)lix p LVj,h{-P) ■ (26d) 


Later we will use similar expressions for contractions with 
other Dirac structures. The neutrino mean-field current itself 
contains spinor contractions of this type and can be expressed 
in terms of the density matrices and pair correlations as 


I?. 

u 



(7 ^ P L)ji t hsPij,sh 

+ (7 ^Pl^IhA^ 


+ (7^ P l) ji,hs K ij,sh 

_|_ (7 ^ P l) sh — Pij,sh) 


(27) 


Notice that in this case there is no implied summation over i 
and j. The fourth term contains a divergent vacuum contribu¬ 
tion that must be renormalized. 

We finally work out the spinor contractions explicitly to 
lowest order in neutrino masses. To this end we introduce 


ra M = (l,p) , •nP = ( 1 , — p) , e M = ( 0 ,e), (28) 


{lix p LYA s h 
(7 Mtsh 
( lu. P LYij,sh 
(7 li P L)ij,sh 


C e~^e, 0 


e 2 p 'G fc /i 


2 p 


e 2 p 


_ e +i<t>T±f, I ’ 
C A» e 2 p ""A *) 


0 * 
u e 2p V 


_ e +i0™i 

e 2p 'G 


(29a) 

(29b) 

(29c) 

(29d) 


where we use the notation introduced in Eq. (15). As an ex¬ 
ample, the vv term, Eq. (29a), reads explicitly 



1 

1 

A 

1 

1 

vv = 1 

1 

1 


T—i 

P 

1 

1 


J 


VV 

ij 


r m 1 m 2 m 3 
mi m 2 m 3 
1 mi m 2 m 3 , 


e +i< ^e* . 

A ' 


VV 

ij 


++ 


= 0 . 


'mi mi mi 

m 2 m 2 m 2 
>m 3 m 3 m 3 , 


e^e 


At ) 


These results agree with those obtained in Ref. [44]. 


(30a) 

(30b) 

(30c) 

(30d) 


HI. MAJORANA NEUTRINO 

From a theoretical perspective, it is quite natural for neu¬ 
trino masses to be of Majorana type. In this case, the two he- 
licity states of a given family coincide with the v and v states, 
the mass term violates lepton number, and there are no ster¬ 
ile degrees of freedom. We work out the modifications of the 
results of the previous section for the Majorana case, concen¬ 
trating again on technical issues. 


A. Two-point correlators and kinetic equations 

In the Majorana case, the momentum decomposition of the 
neutrino field looks the same as for the Dirac case Eq. (9). 
However, because there are no independent antiparticle de¬ 
grees of freedom, the field mode p has the simpler form 

v%, s {^P) = ai, s (t,p)ui, s {p) + uj s (f, -p)vj, s (-p) • (31) 

The creation and annihilation operators satisfy the same anti¬ 
commutation relations Eq. (11) and the bispinors are the same 
as in the Dirac case. 
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The definitions of the two-point correlation functions are 
different because of the different particle content. 


as in Eq. (14). The components of the effective Hamiltonian 
now read [44] 


(2n) 3 S(p-k)p lj!Sh (p) = (a] ift (+k)a i;S (+p)) , (32a) 

(2n) 3 S(p-k)p ijtSh (p) = (at s (-p)oj |h (-k)) , (32b) 

(27r) 3 5(p-k)K ijiS/t (p) = (a j ,/ l (-k)a i)S (+p)) , (32c) 

(2tt) 3 5(p— k)4j,sh(P) = <aj, h (+p)at.(- k )> , (32d) 


KU (p) = r Zh (p) - r f tha (-P ), (38a) 

Klh (p) = r tsh (p) - r t hs (-p) , 08b) 

h^(p) = r% sh ( P ) - rf^ s (- P ) , (38c) 

H^(p) = r?? >fc (p) - r£ fca (-p) . (38d) 


where all operators are taken at the same time t. In the Dirac 
case, Z has no additional information relative to k. Here we 
have additional redundancies 


Pij,sh(ti P) — Pji,hs(ti P) i (33a) 

P) P) > (33b) 

which reflect that Majorana neutrinos have half as many de¬ 
grees of freedom as Dirac ones. Note that in the Majorana 
case, the pair correlations violate total lepton number. 

The mean-field Hamiltonian, bilinear in the neutrino cre¬ 
ation and annihilation operators, has the same form Eq. (7) as 
in the Dirac case. However, as we will demonstrate below, the 
kernel has a more general structure, 

r« = 7 >jPlV f 3 + P L V* + P R V$ . (34) 

The first piece, V- 1 -, is defined as in Eq. (24). In addition, there 
are two scalar pieces 

V* ,R = V2G F l(f , (35) 

depending, as we will see, on the left-chiral and right-chiral 
neutrino mean-field scalar background 

I^ R = {*iPL,RVi) • (36) 

These scalar pieces are missing in the previous literature. 1 
Their explicit form in terms of the density matrices and pair 
correlators will be given in Eq. (41). 

The mean-field Hamiltonian can be written in a form simi¬ 
lar to Eq. (13), 


flmf = 


' p ,s,h L 


0 M(P) r i' J >h(P) a j,h(P) 

+<>) r rjU(p)4fc(-p) 

+ aiA-P) T ij,sh(p) a iAP) 

+ ai,s(-P) T ij,sh(p) a ],h(-P) 


(37) 


where the matrices T l,v , T l,v , etc. are the spinor contractions 
defined in Eq. (12). Using the Heisenberg equation of motion 
with the Hamiltonian Eq. (37) one recovers the equation of 
motion i R = [H, R], where R and H have the same structure 


1 In a private communication, the authors of Ref. [44] agree that these terms 
should indeed be present in the Majorana case. Of course, the presence of 
these terms does not modify the overall structure of the kinetic equations. 


Not all of these components are independent. In particular 

H^(P) = -H^ 5 (-P), (39a) 

H^(P)=-H^%(-P), (39b) 

so only two of the four submatrices of H are independent. 


B. Neutrino-neutrino mean-field Hamiltonian 

The Majorana neutrino interaction with matter is described 
by the same charged- and neutral-current Hamiltonian den¬ 
sities Eqs. (16) and (19) which lead to the same mean-field 
currents of electrons and nucleons—see Eqs. (17) and (20). 

The neutrino-neutrino interaction in the Majorana case is 
also described by Eq. (21). However, Majorana neutrinos vi¬ 
olate lepton-number conservation, and in addition to the four 
lepton-number-conserving combinations considered in Sec. II 
one should also take into account the lepton-number-violating 
combinations (yiVj) and (pMj) which were not included in 
the previous literature. 

To calculate these additional contractions, we use the defi¬ 
nition of the charge-conjugate field v c = C P T , where C is the 
charge-conjugation matrix which has the property C T C = 1. 
Using this definition v = ( v c ) T C and u = C 7 °(i/ c )*, which 
further implies v^P^v = —v c j i1 Pri' c . Therefore, we 
can rewrite the Hamiltonian as — Pl^i] [PjlijPR^j] in 
Eq. (21). TheFierzidentity [55] {^P l )[^Pr]= 2 (Pr][Pl) 
further allows us to rewrite it as 2 [z/iP/jz/j] [zPP^m,] , where 
another sign change was induced by anticommuting the neu¬ 
trino fields. Taking the expectation value of one of the square 
brackets we obtain for the new contribution to the mean-field 
Hamiltonian density 

HZ - V2G f J2 ( [»iPRvf\ lij + \P?PLVi\tf) ■ (40) 

ij 


These new terms supplement the expression for the effective 
Majorana Hamiltonian obtained in the previous literature [44]. 
In Appendix B we reproduce this result using two-component 
notation. 

Two comments are in order here. First, for Majorana fer¬ 
mions v c = v and therefore the resulting contribution to the 
kernel reduces to the last two terms in Eq. (34), while the defi¬ 
nition of left- and right-chiral neutrino backgrounds reduces to 
Eq. (35). Second, DPrv c = and ^P^v = where 
vl = I’[V, which are nothing but components of the Majo¬ 
rana mass term. 
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C. Components of the Hamiltonian matrix H 


The new contributions stemming from neutrino-neutrino in¬ 
teractions can be expressed in terms of the (anti)particle den¬ 
sities and pair correlators. 



{Pl) ji,hsPij,sh + (Pl) ji hs^ij^sh Pij,sh ) 


(PL)ji,hs K ij,l is + (PL)ji,hs K tj,sh 


(41) 


where we have again suppressed the common arguments p 
and (f, p). The notation for the scalar contractions (Pl)™ sh 
etc. is analogous to Eq. (26), except that now there is no "/ ll 
included. 

To lowest order in the small neutrino masses we find, using 
the explicit form of the chiral spinors of Appendix A, 


{Pl)ZsH 
(p T y? 

V-* Ljfij^sh 

(Pl)%sH 
( Pl)Z h 



(42a) 

(42b) 

(42c) 

(42d) 


The components of (Pr) can be obtained from these results 
using the relations (P R )^ sh = [(P L )% h J* and (Pjt)^ ah = 
[(Pl)J7/j S ]*> as we H as similar relations for the remaining two 
components. 

Using the definitions Eq. (38) combined with Eq. (29) and 
the corresponding definition for the scalar case we obtain for 
the vv component of H 


that are linear in the neutrino masses, and additionally terms 
quadratic in the neutrino masses which we neglect here. 

A peculiar feature of Eq. (44) is that its last two lines con¬ 
tain terms proportional to k and re' that are not suppressed 
by the neutrino masses and therefore do not vanish when we 
set the masses to zero. This is somewhat surprising because 
we expect that Dirac and Majorana neutrinos are equivalent 
for m„ —> 0. Therefore, the components of H must coin¬ 
cide in this limit. We return to this question in Sec. IV, where 
we study the case of massless two-component neutrinos and 
demonstrate that in the massless limit these additional terms, 
which are proportional to the lepton-number-violating corre¬ 
lators, are not produced if they are zero initially. 

On the other hand, one important finding of our paper is that 
for a Majorana neutrino with an arbitrary small mass, lepton- 
number-violating correlators are automatically produced and, 
in turn, induce the additional scalar background terms of the 
mean-field Hamiltonian which then affect the dynamics of the 
density matrices. 


IV. WEYL NEUTRINO 


In the previous section we have found that the additional 
scalar contributions to the mean-field Hamiltonian, that natu¬ 
rally arise for Majorana neutrinos, do not vanish in the mass¬ 
less limit. This is somewhat surprising because we expect no 
difference between Dirac and Majorana neutrinos in this case. 
To clarify this paradox we study a single generation of mass¬ 
less neutrinos. The equations presented in this section will 
also be used later to study particle-antiparticle coherence. 


A. Standard two-point correlators and kinetic equations 


H?£h(p) = Seh&ijEi 

+ (7„^&(p)^ - Mfihsi-PKi 
+ (Pl)%sh(v)V* - (Pl)%h.(-p)V£ 

+ (Pn)%sh(PWj ~ (PR)ii,hs (-P)^i • (43) 


In the Weyl case, the momentum decomposition of the neu¬ 
trino field looks the same as for the Dirac case Eq. (9). How¬ 
ever, because a Weyl fermion has only two degrees of freedom 
the field mode p does not carry a spin index, 

P) = a(t,p)u_(p)+b' f (t,-p)v + (-p ). (45) 


The second line generalizes the Dirac result of Eq. (25a) to the 
Majorana case and has been obtained in Ref. [44]. The third 
and fourth lines stem from the the contractions Eq. (40) and 
supplement the previous results. 

The vv term follows from the identity Eq. (39). For the vv 
component we find 

hSJp) = in,PL)Zship)v^ - (7^)£ ha (-P)v£ 

+ (Pl)Zh(pW? - i p L)ths{-P)Vg 
+ {PR)ilsh{v)Vi L i - (Pr)j?M-P) y j L i • (44) 

The vv component follows from replacing vv with vv every¬ 
where in this result. 

An inspection of Eqs. (41) and (42) shows that the last 
two lines of Eq. (43) contain terms proportional to re and re' 


It is automatically left-chiral because the right-chiral compo¬ 
nents of the chiral spinors u_ (p) and v + (—p) vanish in the 
massless limit, see Appendix A. 

If we require lepton-number conservation then the only cor¬ 
relators that we can define are 

(2tt) 3 (5(p - k)p--(p) = (a t (k)a(p)}, (46a) 

(2tt) 3 5(p - k)p ++ (p) = <fo t (—k)6(—p)), (46b) 

(2tt) 3 J(p - k)/c. + (p) = (&(—k)a(p)), (46c) 

(2tt) 3 J(p - k)4-(p) = (« t (p)& t (-k)) ■ (46d) 

Note that we keep helicity indices in these definitions to dis¬ 
tinguish the lepton-number-conserving correlators from the 
lepton-number-violating ones, which we introduce below. 
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We can extract the explicit form of the kinetic equations for 
these correlators from Eq. (2), 

ip- = H Zp- - p-MZ + HZ 4 - - k-+ Hi- , (47a) 

ip ++ = H Zp ++ - - P++ HZ - H Zk. + + 4- W + p , (47b) 
ik. + = HZK-+ - «. + H Z - H^p ++ - p- 

+ H_4, (47c) 

where we omit the arguments (i, p), which are common to all 
the functions, to shorten the notation. Note that for a single 
neutrino generation the first two terms in Eqs. (47a) and (47b) 
cancel each other and we have retained them only to keep the 
resemblance with the general form of the kinetic equations. 

A peculiar feature of Eq. (47c) is that re, i.e. the coher¬ 
ence between 100) and 111) states, is automatically induced 
provided that the mean-field Hamiltonian H has nonzero off- 
diagonals. The off-diagonals can be induced even if all neu¬ 
trino two-point functions are zero initially by, for instance, a 
transverse neutron current. 

The explicit form of the mean-field Hamiltonian can be ob¬ 
tained from Eq. (29) by setting the masses to zero. 


H Z (p) = E + V° - pV , (48a) 

H^(p)=— e*V, (48b) 

H4(p) = —eV, (48c) 

H£! (p) = -E + y° + pV , (48d) 


where E = |p|. Note that the pV term in Eq. (48a) accounts 
for the enhancement (suppression) of the mean-field potential 
for the matter flowing antiparallel (parallel) to the neutrino 
momentum. This has been pointed out in Ref. [34]. 

It remains to express the neutrino current I 1 ' in terms of the 
density matrices and pair correlations. For its time component 
we obtain from Eq. (27), 1° = f p i, where £(t, p) = p(t, p) — 
p(t, —p) has the meaning of lepton number in mode p. For 
the spatial components we find 



which coincides with the result of Ref. [44]. 

B. Lepton-number-violating correlators and kinetic equations 

If we allow for (////) and (vv) contractions then, similarly to 
the Majorana case, the mean-field Hamiltonian receives con¬ 
tributions of the type Eq. (40). Because Weyl fields satisfy the 
condition P^is = v we can rewrite Eq. (40) as 

HZ = V2G f J2 (M I L + M I R ) (50) 

(see Sec. Ill and Appendix B for more details) where now 

I L = QFv) and I R = (vv c ). (51) 

As has been mentioned above v r v and uu c have the struc¬ 
ture of the Majorana mass term, which is known to violate 


lepton number. Therefore, we expect that also for the Weyl 
neutrino the mean-field Hamiltonian Eq. (50) leads to lepton 
number violation. However, for Weyl neutrinos the inclusion 
of these additional terms is somewhat artificial because, as we 
show below, these correlations are not produced if they are 
zero initially. They are considered here to better understand 
the Majorana case, where they are naturally produced by the 
lepton-number-violating interactions. 

The contribution of Eq. (50) to the mean-field Hamiltonian 
is given by 

H mf = J [ a t( p )rrf(p)at(-p) 

+ & t (p)r^(p)6 t (-p) 

+ a(-p)n / _ ! '(p)a(p) 

+ 6(-p)r^(p)6(p)], (52) 

and strongly resembles the mean-field Hamiltonian of Majo¬ 
rana neutrinos Eq. (37). From the structure of Eq. (52) it is 
evident that, as expected, it leads to the violation of lepton 
number. To take this into account we are forced to introduce 
the following lepton-number-violating correlators, 

(2tt) 3 <5(p - k)K-(p) = (a(-k)a(p)) , (53a) 

(2tt) 3 <5(p - k)K ++ (p) = <6(—k)6(p)> , (53b) 

(2tt) 3 5(p - k)4_(p) = (a t (p)a t (—k)) , (53c) 

(2tt) 3 (5(p - k)4+(p) = {b\p)b\-k)} , (53d) 

which also resemble the Majorana definitions Eq. (32). These 
correlators are dictated by the structure of the Hamiltonian 
Eq. (52) and are the only lepton-number-violating correlators 
we consider in this section. If we wanted to consider all other 
possible correlators we would be back to the Majorana case 
with zero neutrino masses. 

The lepton-number-violating correlators contribute to the 
dynamics of the lepton-number-conserving ones, 

ip- = ... 

+ (TZ - [VZ] T ) kL - K— (TZ - [TZf) , (54a) 

ip++ = ■■■ 

- (r£ - [4:] T ) K++ + 4 + (r "Z - [r^H, (54b) 

where ellipses denote terms on the right-hand side of Eq. (47), 
and the superscript T stands for transposition of the flavor and 
helicity indices, as well as inversion of the momentum. Com¬ 
paring Eq. (54) with Eqs. (38b) and (38c) we see that we auto¬ 
matically recover the “Majorana” definitions of the Hamilto¬ 
nian matrix. Note that to avoid confusion with the definitions 
of the elements of the mean-field Hamiltonian, which are dif¬ 
ferent for Dirac and Majorana neutrinos, we write the right- 
hand side of Eq. (54) directly in terms of spinor contractions 
defined in Eq. (12). The dynamics of n.+, see Eq. (47c), does 
not receive any corrections. The reason is that the components 
of the mean-field Hamiltonian needed to form the right spin 
combination with the lepton-number-violating correlators in 
Eq. (47c) are zero for Weyl neutrinos. The kinetic equations 
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for the lepton-number-violating pair correlations read 

ik- = T™k- - K-(-[TT] t ) - p-(Tl 9 - [TT] t ) 

- (rf - [Pf] T ) [ P ..] T + (r pp - [rr_T) , (55a) 

ik ++ = (-[r pp ] T )« ++ - K++ r pp - {TZ - KI] t )p+ + 

- [p ++ ] T (vii - inif) + in: - [r::n . (55b } 

Their form can be guessed from Eq. (47c) by replacing com¬ 
ponents of the mean-field Hamiltonian with their “Majorana” 
counterparts, taking into account that ITT = I'H = 0, and re¬ 
placing p- by [p-] T as well as p ++ by [p++] T ■ 

Using the explicit form of the chiral spinors (see Appendix 
A) we obtain 


TlH (p) = +e +i( ^V L , 

(56a) 

Tll(p) =-e~^V R , 

(56b) 

r_T(p) = 

(56c) 

r p [[(p) = —e +i ^V L , 

(56d) 

where V ,A,l] = \f2GyI L(H> are defined analogously to 
Eq. (35). Let us now recall that I L and I R are produced 
only by neutrino self-interactions and are proportional to the 
lepton-number-violating pair correlations. 

I L = [ e-^[«_ , 

Jp 

(57) 


and a similar expression for I R . Thus, if the lepton-number- 
violating correlators are zero initially, then the components in 
Eq. (56) are zero and k~ and k++ remain zero in the course 
of the system’s evolution. For this reason for Weyl neutrinos 
the inclusion of lepton-number-violating correlators is rather 
artificial because they could only exist if they were put in by 
hand initially. 

This observation explains why similar contributions do not 
vanish for Majorana neutrinos in the limit of zero neutrino 
masses. While such lepton-number-violating correlators can 
be introduced by hand as an initial condition, they can dynam¬ 
ically evolve only in the presence of a nonvanishing Majorana 
mass. 


V. ELECTROMAGNETIC BACKGROUND FIELDS 

A supernova environment is characterized not only by mat¬ 
ter currents, but also by strong magnetic fields. Electromag¬ 
netic fields polarize both background media and the vacuum. 
Although neutrinos do not couple directly to the electromag¬ 
netic fields, they feel the induced polarization. The coupling 
to a polarized background medium has been treated in the pre¬ 
vious sections. We now turn to the interaction with the vac¬ 
uum polarization. 

The effect of vacuum polarization is described by electro¬ 
magnetic form factors. The most prominent examples, the 
magnetic and electric dipole moments, are inevitable for mas¬ 
sive neutrinos and have to be included to obtain consistent 


evolution equations linear in the neutrino mass. The main ef¬ 
fects of electromagnetic fields are spin and spin-flavor oscilla¬ 
tions, which can be significant. We treat Dirac and Majorana 
neutrinos separately. 


A. General vertex structure 

The coupling of neutrinos to an external vector potential A 1 ' 
can be written as an effective vertex TL em = A^vT^v, where 
T^ contains all irreducible combinations of Lorentz vectors 
and pseudovectors generated by external momenta and Dirac 
matrices. Neglecting a hypothetical minicharge, in coordinate 
space the most general Hamiltonian density can be reduced to 

+ d v F^ u Vi(f% 7 m 75 )vj , (58) 

where the electromagnetic field-strength tensor is defined as 
usual, F I1V = d fl A v - and a' ,v = f [7T7"]. The 

form factors are Jm (magnetic), Je (electric), /q (reduced 
charge [28]), and Ja (anapole). The form factors carry gen¬ 
eration indices. Diagonal elements describe the usual elec¬ 
tromagnetic properties of a neutrino in the mass basis, and 
reduce to electromagnetic moments in the static limit. The 
off-diagonal elements describe transitions between neutrinos 
of different masses. Some components of the Hamiltonian 
matrix have been calculated in Refs. [28, 56]. 

Maxwell’s equations tell us that d v F llv = —J R m , where 
Jem is some charged matter background that sources electro¬ 
magnetic fields. In supernovae, the sources are electrons and 
protons. In the Standard Model with massless neutrinos, the 
value for the anapole moment has to be f a = —/q to repro¬ 
duce the left-chiral form of the interaction. For models with 
neutrino masses, the Hamiltonian matrix might obtain contri¬ 
butions that are not purely left-chiral, but we assume that these 
are always small so that we can neglect them. The charge 
and anapole form factors then only yield radiative corrections 
to the left-chiral tree-level coupling in Eq. (16). We neglect 
these moments because we are not interested in corrections to 
leading-order effects. However, for completeness, we give the 
spinor contractions for right-chiral currents in Appendix C. 


B. Dipole moments of Dirac neutrinos 

To study the dipole moments, we first turn to the some¬ 
what simpler case of Dirac neutrinos. A Dirac neutrino has 
diagonal magnetic and electric moments. Because we assume 
neutrinos to carry no charge, // = /m(0) is defined as the 
magnetic moment and e = /b( 0) as the electric dipole mo¬ 
ment [28]. In the minimal extension of the Standard Model, 
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the magnetic moments are found to be [57] 


Fij ~ 


3 eV2G F (rrii + rrij) 

2(4^p v 

. 3e\/2GF , . ( m 2 




■Fij 1 


E u l(^) U °1’ 


a—e,fi,T 


(59a) 

(59b) 

(59c) 


where m T is the tau mass. Note that the electric dipole mo¬ 
ment does not have a diagonal component because it would 
violate CP [28], and that the transition electric dipole moment 
carries a phase relative to the transition magnetic dipole mo¬ 
ment. Numerically, the above expressions yield for the diago¬ 
nal magnetic moments 

IHi - 3.2 x lO" 19 /i B , (60) 

where // B is the Bohr magneton. The transition moments are 
Hii - -3.9 x 10- 23 J-y ; (6 i a ) 

ey ~ 3.9i x 10- 23 J-y Mb ■ (61b) 

Note that the transition moments are much smaller than the di¬ 
agonal moments due to Glashow-Iliopoulos-Maiani suppres¬ 
sion. 


C. Hamiltonian matrix for Dirac neutrinos 


where we identify H vv = T vv , W vv = F vv , etc., and the mi¬ 
nus sign in the metric g^ v = diag (1, — 1, — 1, —1) has already 
been taken care of. In Eq. (62), the form factors still depend 
on the momentum transfer. For the vv and vv components, 
the form factors contain q^ = p^ ut — pM, where —> 0 in the 
forward-scattering limit. These components are then propor¬ 
tional to the dipole moments. For the neutrino-antineutrino 
components of the H matrices, the argument of the form fac¬ 
tor contains l 2 with l M = pg Ut + p l ‘ n , the sum of neutrino and 
antineutrino momenta. In the forward-scattering limit this re¬ 
duces to l 2 = (2 E) 2 , and the dependence of the form factors 
on the four-momentum is important. 

The coupling of the magnetic field to the electric form fac¬ 
tor (superscript eB) is 


HSX=-(iV > 7)^/g(« 3 )B, (63a) 

HSt=-(i7°7)J, k /2(i 2 )B, (63b) 

Hgt=-((7°7)" k /!(! 2 )B, (63c) 

HgS=-(>7°7)" k /g(9 2 )B. (63d) 


The coupling of an electric field to the magnetic form factor is 


Sr FhFiZjtitf )e, 

(64a) 

St =(‘7°7 

(64b) 

Sr =(‘7 0 7)“ „S(! 2 )E, 

(64c) 

St =C°7 

(64d) 


We treat electromagnetic effects on the same footing as 
background matter. To this end, we have to evaluate the com¬ 
ponents of the Hamiltonian matrix, which, for Dirac neutri¬ 
nos, are equal to the spinor contractions in Eq. (12). For the 
contractions, we need to evaluate the Lorentz structure of the 
vertex in Eq. (58). 

Considering only magnetic and electric form factors, 
the Hamiltonian reduces to 75 ) 27 -, 

which depends on the electric and magnetic fields, E and B, 
through F ,,v . The Lorentz structure can be decomposed into 
the contractions (* 7 ° 7 ) ; sh and ( 7 ° 775 ) i ■ h , the latter ap¬ 
pearing through the identity e abCr y 0 'y Cr y 5 = a ab with spatial 
indices a,b,c = 1,2 or 3, and the asymmetric tensor e abc . 
These contractions are three-vectors that are contracted with 
the electric and magnetic fields. We calculate the contractions 
in momentum space. 

Explicitly, the coupling of the magnetic field through the 
magnetic form factor (superscript // 11 ) has the structures 

H St = - ( 7 ° 77 s)) 7 l /jJ(<! 2 )B, (62a) 

HSS = - (7°77a)J, 1 /SC 2 )B, (62b) 

K’S = - , (620 

HS“‘’=-(7°77 5 )"7S(9 2 )B, (62d) 


which is indicated by p, E, and to the electric form factor, eE, 

H^=-(7°775);^/|(9 2 )E, (65a) 

H Zl = - (7°775);. sh f£0 2 )^ (65b) 

H^=-(7°775)^/|(( 2 )E, (65c) 

H^=-(7°775);^/|(9 2 )E. (65 d) 

One can see that a magnetic field couples to both, the elec¬ 
tric and the magnetic form factor. Also electric fields cou¬ 
ple to both form factors. This can be understood as follows. 
In the neutrino rest frame, the magnetic field only couples to 
the magnetic dipole moment, and the electric field only cou¬ 
ples to the electric dipole moment (if any), as suggested by 
the nomenclature. Lorentz covariance then demands that both 
electric and magnetic fields couple to the magnetic form fac¬ 
tor in a system where the neutrino moves with nonzero veloc¬ 
ity. A moving neutrino also exhibits spin precession in a pure 
electric field through its magnetic moment [58]. 

The Lorentz structure of Eqs. (62)—(65) can now be readily 
calculated. In contrast to the previous sections, we neglect all 
contributions proportional to the mass since the magnetic and 
electric form factors are small and, in the models considered 
here, proportional to the neutrino mass already. The ( 7 ° 775 ) 
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components are 


(7°775) 

(7°775) 

(7°775) 

(7°775) 


VV ^ 

( 0 e 

-H0g* \ 

ij,sh 

e-^e 

0 ) ■ 

vv 

( -e +i ^p 

0 

ij,sh 

l o 

~e~^p 

vv ^ 

<Oh 

T 

0) 

1 

0 

ij,sh 

V 0 

^g+*0p 

vv ^ 

( 0 


ij,sh 

l —e+^e 

0 


(66a) 

(66b) 

(66c) 

(66d) 


The remaining Lorentz structures are of the form ('<7°7). 
They read 


(*r°7)£. fc 

(*7°7 )t h 


0 ie+^e* \ 

-ie-i't’e 0 ) ’ 

p 0 \ 

0 ie-^p ) ’ 


(*r°7)S:. fc 

(*7°7)" h 


ie 0 \ 

0 -ie +i ^p ) ’ 

0 ie-^e* \ 

-ie+^e 0 J ‘ 


(67a) 

(67b) 

(67c) 

(67d) 


neglecting the model-dependent amplitudes, one can compare 
the moments of Dirac and Majorana neutrinos. The main dif¬ 
ferences are that the Majorana amplitudes contain Majorana 
PMNS matrices, which may contain more phases than Dirac 
PMNS matrices, and that Eq. (68) has to be taken into account 
for Majorana neutrinos. 

After these adjustments, the off-diagonal form factors of 
Majorana neutrinos can be obtained from Eq. (59a). One finds 
that they depend on the relative CP-phases of two neutrino 
species [59], The relative phase can either be equal or oppo¬ 
site, i.e., the ratio is ±1. For neutrinos with equal CP-phases, 
the magnetic transition moments vanish [59], while for oppo¬ 
site CP-phase the magnetic transition moments are nonzero 
and can be obtained from Eq. (59a) by substituting Tj with 
c Xi\mTij. 

For electric dipole moments, the role of the CP-phases is 
inverted. Opposite CP-phases force the electric transition mo¬ 
ments to vanish, while for equal CP-phases the electric tran¬ 
sition moments are nonzero and are obtained by substituting 
with 2R eTij [59] in Eq. (59b). 


E. Hamiltonian matrix for Majorana neutrinos 


To this level of approximation, the vv and vv components 
are diagonal in helicity space, i.e., electric and magnetic fields 
mainly couple spin-0 neutrino-antineutrino pairs. Because the 
diagonal is proportional to p, the relevant field components 
are those parallel to the momentum of the neutrinos. The vv 
and vv components are off-diagonal in helicity space. The 
dominant effect of magnetic and electric fields on neutrinos 
and antineutrinos is spin precession. Here the transverse com¬ 
ponents of the electromagnetic fields contribute. The longitu¬ 
dinal components enter on the diagonals in the next order of 
the expansion in m/E and are therefore omitted. 


The density matrix formalism naturally reproduces the re¬ 
sults for the electromagnetic moments discussed in the last 
section. Similarly to Eq. (68), each component of the Hamil¬ 
tonian matrix has two contributions from T contractions, e.g., 

H ij>fc(p) = r i J -,ah(p)- I 'ji,hs(-p ^ see E q- ( 38 >- The S P in0r 
contractions r i/l/ and I ' v,/ have the same structure as for Dirac 
neutrinos; see Eq. (12). Again neglecting the model depen¬ 
dence, the only difference is that the Dirac PMNS matrix have 
to be replaced by the Majorana PMNS matrix. For example, 
a magnetic field coupling to a Majorana neutrino through the 
magnetic form factor yields 


D. Dipole moments of Majorana neutrinos 

For Majorana neutrinos, electromagnetic transitions always 
contain two contributions, e.g., 

(''pj^ I^Pin) = (^PoutE^Pin _ ^PinT^fp^,) • (68) 

This difference of two amplitudes leads to the cancellation of 
all the diagonal moments except for the anapole moment [28]. 
This can also be understood by noting that the last two terms 
(including the minus sign) in Eq. (68) are charge conjugates 
of each other. Because the Lorentz structure of the magnetic, 
electric, and charge form factors are C-odd the combination 
vanishes. The Lorentz structure of the anapole moment is el¬ 
even and does not cancel. 

Because the magnetic moment of the Majorana neutrino 
vanishes, it does not couple directly to a magnetic field. How¬ 
ever, magnetic fields polarize the background medium, and 
this effect does lead to helicity oscillations, see Sec. III. 

Electromagnetic moments of neutrinos depend on the de¬ 
tails of the mechanism that creates the neutrino mass. When 


H^. + = -[/ii(^)-c.c.)]e + ^B 
= - 2iIm[/^(g 2 )]e+^e*B , 

where we have used the Hermiticity of the form factors. In 
the static limit, 2/'Tm[/^] is the magnetic moment of Majo¬ 
rana neutrinos [57]. It is zero for equal CP-phases since Tij 
becomes real. It is nonvanishing for opposite CP-phases be¬ 
cause Tij becomes imaginary. An analogous argument holds 
for the electric dipole moment. 


VI. HELICITY COHERENCE 

In this section we neglect pair correlations and discuss he¬ 
licity coherence effects. To separate the latter from the usual 
flavor coherence effects, we consider only one neutrino gener¬ 
ation. Furthermore, for definiteness we assume that neutrinos 
are Dirac particles. 
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A. Order-of-magnitude estimate 


of the Hamiltonian matrix. 


Two different mean-field backgrounds cause spin oscilla¬ 
tions and create spin coherence: matter and neutrino currents, 
and electromagnetic fields. However, it is not clear which of 
these is dominant in a supernova. In the following we perform 
a crude estimate. 

For Dirac neutrinos without pair correlations, the kinetic 
equations of neutrinos and antineutrinos decouple, ip = 
[H"", p\ and ip = [H" 1 ', p], and, for one family, we only have 
to look at a 2 x 2 subsystem of the full evolution equation. 
We start with a matter background with nonrelativistic veloc¬ 
ity j3, which flows orthogonal to the neutrino’s momentum. 
The Hamiltonian matrix reads 

( 1 — b\ 

^ V (§/? 0 J ’ ™ 

where V is the usual matter potential. For instance for 
or u T it is given by V = G-pn n l \/2, where n n is the neutron 
density. We have omitted the neutrino kinetic energy because 
it is diagonal in helicity space, and for a single generation 
trivially cancels in the commutator. The dependence of the 
diagonal terms of the Hamiltonian on the parallel flux and the 
dependence of the off-diagonal terms on the orthogonal flux 
were discussed in Refs. [35, 60]. 

In a derivation similar to the one that leads to Eq. (70), we 
obtain the 2x2 subsystem of the Hamiltonian matrix for a 
neutrino in a transverse magnetic field 

H w » ~I*B J) (71) 

(see Sec. V for more details). Spin coherence is instigated by 
the off-diagonals of Eqs. (70) and (71), and to find the relative 
importance of the matter and magnetic contributions it is suf¬ 
ficient to estimate their relative size. Typical magnetic fields 
in a supernova are of order 10 12 G and much larger in magne- 
tars. Using the standard value for the magnetic moment given 
in Eq. (60), and assuming a neutrino mass of 0.1 eV, we find 
for the contribution of the magnetic field ///i ~ 10~ 16 eV. 
For a typical neutron mass density 10 12 g/cm 3 , which cor¬ 
responds to a number density n n ~ 10 4 MeV 3 , the matter 
potential is of the order of the neutrino mass, V ~ 0.1 eV. 
Thus, for a typical momentum p ~ 30 MeV we obtain 
Vf3m/(2p) ~ 10 _10 ^eV. For maximal background veloc¬ 
ities of 3000 km/s, /3 ~ 0.01, the matter contribution dom¬ 
inates. Surprisingly, the magnetic field is only important if 
the background moves very slowly, if the matter density has 
decreased sufficiently, or if the magnetic moment is enhanced. 

Turning now to the density matrix, the size of the off- 
diagonal elements depends on the initial conditions and his¬ 
tory of the evolution. To obtain a rough estimate, we can as¬ 
sume that the system has reached equilibrium and, hence, its 
previous evolution is irrelevant. In equilibrium, the system is 
in an eigenstate of the Hamiltonian, i.e., H !/,/ and p commute. 
This condition alone allows us to express the off-diagonals of 
the density matrix in terms of the diagonals and components 


P-+ = 


(n __ P++) 


(72) 


Keeping only the (dominant) matter contribution, Eq. (70), we 
find p-+ = (p- — p ++ )mf3/2p ~ m/3/2p. For m ~ 0.1 eV 
and a typical momentum p ~ 30 MeV this results in p.+ ~ 
10 -11 , where we have used /? ~ 0.01. 

The same result can be obtained by noting that if p and H vv 
commute, they can be simultaneously diagonalized by a rota¬ 
tion that mixes positive- and negative-helicity states. The ro¬ 
tation angle is tan 2$ = m/3/p. Considering e.g. the p— = 1 
eigenstate of the diagonalized Hamiltonian and rotating back 
to the basis where the Hamiltonian has the form Eq. (70) we 
find to leading order 



(73) 


The corrections to the diagonals are not included in Eq. (73) 
because they are of the order of 6p- ~ Sp++ ~ p 2 + ~ 10 -22 
and are therefore negligibly small. 


B. Resonant enhancement 


For a magnetic field, the diagonal elements of the Hamilto¬ 
nian matrix, Eq. (71), are zero for very relativistic neutrinos. 
This allows for the magnetic fields to completely flip the spin 
of a population of neutrinos. On the other hand, the diagonals 
of Eq. (70) are in general nonzero and suppress a complete 
conversion. In general, the matter contribution is given by 


H w 


fv°-V, §v±) 

[%V ± 0 ) 


(74) 


see Eq. (43), where V u = pV and V± = eV are compo¬ 
nents of the matter flux parallel and orthogonal to the neu¬ 
trino momentum. Thus, if there are relativistic currents par¬ 
allel to the momentum of the neutrino such that the diago¬ 
nals vanish, Eq. (72) implies that a resonant enhancement of 
the spin conversion is possible. The possibility to generate 
the spin conversion by an orthogonal flux of matter, and the 
cancellation of the matter effect for relativistic matter moving 
along the direction of the neutrino momentum were first dis¬ 
cussed in Refs. [35, 60] on the basis of the Lorentz-covariant 
Bergmann-Michel-Telegdi equation. In Refs. [36, 38] these 
effects have also been studied using the formalism of nonequi¬ 
librium quantum field theory. In the context of resonant lep- 
togenesis the formation of flavor and helicity correlations in 
medium and the derivation of flavor-covariant transport equa¬ 
tions able to account for helicity correlations has been dis¬ 
cussed in Ref. [61]. 

For vanishing diagonals, Eq. (74) can be rotated into its di¬ 
agonal form with a rotation angle i? = 7t/4. In other words, 
mixing of the helicity states becomes maximal, similarly to 
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the Mikheyev-Smirnov-Wolfenstein resonance mixing, and 
hence in equilibrium 






(75) 


where we have again assumed that the system is in an eigen¬ 
state of the diagonalized Hamiltonian. Outside of the core, a 
supernova is far from equilibrium, but nonlinear feedback can 
enhance the spin-flipping processes [38]. 

Making use of Eq. (24), we can rewrite the resonance con¬ 
dition, - H™ = V° - V H = 0, in the form [38] 



(76) 


where Y e = n e /riB and Y u = (n v — rip)/ns are the electron 
and neutrino asymmetry fractions respectively and ns is the 
baryon number density. The resonance condition can poten¬ 
tially be fulfilled in or near the proto-neutron star in a core¬ 
collapse supernova, or near the central region of a compact 
object merger; see Ref. [38] and references therein. 


C. Lorentz covariance 

Helicity coherence builds up only if the off-diagonal ele¬ 
ments of the Hamiltonian matrix differ from zero. On the 
other hand, because the off-diagonals are proportional to the 
component of the matter flow orthogonal to the neutrino mo¬ 
mentum, one can always find a frame where the ol'f-diagonals 
vanish and no helicity coherence builds up. In other words, at 
first sight physical results seem to depend on the frame. This 
raises the question of Lorentz covariance of the kinetic equa¬ 
tions. 

To be specific, let us consider the following simple exam¬ 
ple. We have two identical observers moving with velocity /? 
with respect to each other. In the frame of the first observer, 
the neutrino has momentum p along the z axis and the matter 
is at rest, i.e. V n = V± = 0, 


» V (J °) . (77) 

Thus no helicity coherence builds up. In the frame of the sec¬ 
ond observer which moves with velocity /3 along the x axis 
the Hamiltonian is no longer diagonal, 

<78) 

and we expect helicity coherence to build up. Here 7 is the 
usual Lorentz factor and V and p denote the potential and neu¬ 
trino momentum in the frame of the first observer. 

Do the Hamiltonian matrices Eqs. (77) and (78) lead to dif¬ 
ferent physical results ? The answer is no, but to demonstrate 
this point we need to take into account that a helicity state is 
also not Lorentz invariant. Let the neutrino be in a state of 


definite helicity in the frame of the first observer, e.g. \pz, -), 
where p is the absolute value of the neutrino momentum and z 
is the unit vector along the z axis. The corresponding density 
matrix reads 

P = If®, -)(pz, “I = q) • (79) 


The Hamiltonian matrix Eq. (77) and the density matrix com¬ 
mute and therefore the latter is constant in time. The boost 
to the frame of the second observer transforms \pz, -) into a 
mixed helicity state with momentum q, 

IV 7 ) = c 6>/2 1 q, -) — se/ 2 |q, +), (80) 

where 6 is the angle of Wigner rotation around y with tan 9 = 
—mp/p. Note that the rotation angle vanishes in the limit of 
zero neutrino mass which reflects chirality conservation. The 
density matrix develops off-diagonal elements, 


P= IV’XV’I 


1 /1 + Cg — Sg 

2 l -sg 1 - Cg 


(81) 


The Hamiltonian matrix and the density matrix again com¬ 
mute. In other words, the second observer sees a mixed he¬ 
licity state which, as expected, is also time independent. This 
result reflects Lorentz covariance of the kinetic equations, the 
lesson being that one has to transform the initial conditions 
consistently to obtain covariant results. 

Let us now consider this result from a slightly different 
viewpoint. In each frame, we can diagonalize the effective 
Hamiltonian by performing a Bogolyubov transformation that 
mixes annihilation (creation) operators of the positive- and 
negative-helicity states, a s —> c$a s + s$a- s . In particu¬ 
lar Eq. (78) is diagonalized by a Bogolyubov transformation 
with the angle tan 2$ = m/3/p. This transformation brings 
the density matrix Eq. (81) back to the form Eq. (79). In 
other words, there is a connection between the Lorentz and 
Bogolyubov transformations. In particular, if in every frame 
we diagonalize the Hamiltonian then the transformed density 
matrix remains invariant under the boosts. 

To summarize, as far as helicity coherence is concerned, 
both the Hamiltonian and the density matrix transform un¬ 
der Lorentz boosts in such a way that the kinetic equation is 
Lorentz covariant. We will rely on this result in the discussion 
of particle-antiparticle coherence whose Lorentz transforma¬ 
tion properties are not as evident as for the helicity coherence. 


VII. PARTICLE-ANTIPARTICLE COHERENCE 

In this section we discuss particle-antiparticle coherence. 
In contrast to helicity coherence, which requires nonzero 
neutrino masses, and flavor coherence, which in addition to 
nonzero masses requires the existence of several neutrino gen¬ 
erations, particle-antiparticle coherence arises already for a 
single massless neutrino generation. As has been discussed 
in the previous section, for a massless neutrino the only “nat¬ 
ural” correlators are p —, p ++ and To shorten the notation 
in this section we suppress the spin indices. 
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A. Quantum-mechanical example 

To clarify the meaning of the particle-antiparticle coher¬ 
ence, let us first study in more detail the simple quantum- 
mechanical example briefly discussed in the Introduction. We 
consider a system that can be in a linear combination of one 
of four pure states. These are i) the empty state 100) with¬ 
out particles; ii) the paired state 111) = a)(p)f^(-p)|00), 
which contains a neutrino with momentum p and an antineu¬ 
trino with momentum -p; iii) the one neutrino state 110); and 
iv) the one antineutrino state |01). Note that in all these states 
the antineutrinos stream in the direction opposite to that of 
neutrinos. A general state can be expressed in terms of these 
four states, \ip) = A O o|00) + An|ll) + Ai 0 |10) + A O i|01), 
where the coefficients Aij are time dependent and normalized 
to unity, |A 00 | 2 + |^n| 2 + |-4io| 2 + |-4oi| 2 = 1. 

In analogy to Eq. (13) we write the Hamiltonian in the form 


neutrinos with momentum p whereas p(t, p) describes an¬ 
tineutrinos with momentum —p, the conservation of p — p 
implies that re induces the production of neutrino-antineutrino 
pairs with opposite momentum. 

The kinetic equation for re describes a driven harmonic os¬ 
cillator with frequency H"" — Y\ vv ~ 2 E. Hence re oscillates 
with twice the neutrino energy as expected. 

From the definition re = Aq 0 Ah we see that nonzero parti- 
cle-antiparticle coherence means that the system is not in an 
eigenstate of the unperturbed Hamiltonian, but instead in a 
mixture of the |00) and j 11) states, i.e., in a squeezed state. 
Such states do not have a definite particle number. This obser¬ 
vation clarifies the physical meaning of particle-antiparticle 
coherence. 


B. Order-of-magnitude estimate 


H = a t (p)H w a(p) + a t (p)H" p 6 t (-p) 

+ 6(-p)H P "a(p) - 6 t ( -p)H^6(-p). (82) 

The Schrodinger equation for the coefficients A r) then splits 
into three independent equations, 



id t A w = H m A 10 , (83b) 

id t A 01 = -H DD A m . (83c) 

Thus the evolution of the single-particle states completely de¬ 
couples because a homogeneous background medium cannot 
mix states of different total momentum. On the other hand, 
the 100) and 111) states have zero momenta and therefore can 
be mixed by a homogeneous medium through the H vv term of 
the Hamiltonian. However, the |00) and 111) states have dif¬ 
ferent angular momentum. Hence, an anisotropic background 
medium, e.g. a transverse matter flux, is needed to absorb the 
angular momentum and to mix the two states. 

To make the connection to the density matrix equations, we 
note that the number of neutrinos and antineutrinos is given 
by p = |An| 2 + |Ai 0 | 2 and p = |Au| 2 + |A 0 i| 2 respectively. 
Their time evolution can be derived from Eq. (83) and takes 
the form expected from Eq. (47), 

p = -2Im(H pl 're) , (84a) 

p=- 2Im(hr l 're) , (84b) 

if we identify re = Aq 0 Ah. Equation (83) also leads to an 
evolution equation for re 

ik= (H uv -H DD )K + H vi> {1- p-p) , (85) 

which can be obtained by using the normalization of the state 
| i/j). Equation (85) is again consistent with Eq. (47) and coin¬ 
cides with the result of Ref. [44] in the one-flavor limit. 

From these kinetic equations we can infer that while p and 
p are not separately conserved in the presence of nonzero re, 
their difference is conserved [44]. Because p(t, p) describes 


As a next step we perform an order-of-magnitude estimate 
of re. For a single neutrino generation the extended density 
matrix reduces to a 2x2 matrix of the form, see Eq. (1), 

" = < 86 > 

where now p and p are real numbers and re is a complex num¬ 
ber. We again start with an example of a matter background 
with nonrelativistic velocity /?, which flows orthogonal to the 
neutrino’s momentum. Then, as follows from Eq. (48), the 
Hamiltonian matrix reads, see also Eq. (3), 

h=b (J-°i) +v (-W)- <87) 

Unlike for helicity coherence, the neutrino kinetic energy E 
no longer cancels out in the commutator. 

Similarly to the case of spin coherence we can get a crude 
estimate of the re magnitude by assuming that the system has 
reached equilibrium and hence re = 0. Equation (85) then 
gives 

\\vD 

re = ——- tt=( 1 - p — p). (88) 

If we insert this result into Eq. (84) and use the Hermiticity of 
the Hamiltonian matrix, we see that indeed p = p = 0. 

Let us assume for a moment that the neutrino-neutrino in¬ 
teractions are small compared to the neutrino interaction with 
matter. For typical supernova parameters V ~ 0.1 eV and 
E ~ 30 MeV and we then find V/E ~ 10 -9 . Thus to a good 
approximation - H p *) ~ V/3/2E ~ 10“ n , where 

we have used (3 ~ 0.01. Because typically \1 — p — p\ ~ 1 
we conclude that the “natural” size of the particle-antiparticle 
coherence isre~10 -11 . 

The same result can be obtained by noting that in equi¬ 
librium R and H commute and can be simultaneously diag¬ 
onalized by a Bogolyubov transformation that mixes neutri¬ 
nos of momentum p with antineutrinos of momentum — p. 
Under this transformation the creation and annihilation op¬ 
erators transform as a(p) —> e~ l ^ 2 c& a(p) + e l ^ 2 Si) (—p) 
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and b^(—p) — y e l ^ 2 c$ p) — e l ^ 2 sg a(p) respectively, 
where the phase </> = arg H"" and the rotation angle is given 
by tan2$ = 2\H vS> \/(H vv - H pp ) ~ V/3/E. In the basis 
where the Hamiltonian is diagonal, the system is described by 
(anti)neutrino densities, which we denote by g and g respec¬ 
tively, and a pairing correlator, which we denote by x. From 
the transformation properties of the creation/annihilation op¬ 
erators, we can infer the following relations 


p = c#g - CfiSfiX - + 4(1 - q) , (89a) 

P = c#Q ~ c#s$x - c.sSdH* + 4(1 - Q), (89b) 


K = e** 


4 >t + CtfSdQ - C$S#(1 - 0) ~ 4 ^ 


(89c) 


see Ref. [43] for a detailed discussion. Eigenstates of the diag¬ 
onalized Hamiltonian are characterized by x = 0. Assuming, 
e.g., that the system is in an eigenstate of the diagonalized 
Hamiltonian with some g and g, and rotating back to the ba¬ 
sis where the Hamiltonian has the form Eq. (87), we find to 
leading order 


depends on the pair correlations of all other momentum modes 
q. Note also that the phase-space integral in Eq. (91) is un¬ 
bounded. Pairing correlations with a momentum typical for 
the supernova environment couple to pairing correlations of 
arbitrary high momentum. This pushes us beyond the limi¬ 
tations of the Fermi approximation, and in principle a fully 
renormalizable theory has to be studied to make sense of these 
high-momentum modes. To stay within the realm of applica¬ 
bility of the effective theory, we use a phenomenological cut¬ 
off |q| = Mw in the phase-space integrals. 

To estimate the contribution of the k terms to the integral in 
Eq. (91), we take into account that pair correlators of different 
momentum modes oscillate incoherently such that we can re¬ 
place k by its approximate mean value, k « Y\ vv /2E, where 
we use that V <C E and assume p + p <C 1 in Eq. (88). To 
proceed we recall that W v = — e*V, see Eq. (48), where V 
is the total potential that includes matter and neutrino contri¬ 
butions. Note further that V is momentum independent. With 
these substitutions, the integrals involving k in Eq. (91) read 


R ~ 


( g 

v£ 

\ IE 



10 - n \ 


(90) 



V2G f M 2 v 
2 3tt 2 ’ 


(92) 


which again leads to the tiny k ~ s# ~ 10 -11 . 

Pair correlations themselves are not measurable, and only 
their effect on the number densities can be observed. A quick 
inspection of Eq. (89) shows that in equilibrium the difference 
between e.g. p and g is of the order of 4 ~ k 2 . In other 
words, the induced corrections to p and p are quadratic in k. 

This can also be understood from Eq. (85). If the system 
has not yet reached equilibrium, then k oscillates around its 
stationary value Eq. (88), provided that the components of the 
Hamiltonian matrix only vary slowly with time compared to 
H"' y — This assumption allows us to approximate the 

evolution of n as a driven harmonic oscillator with an am¬ 
plitude that depends on the initial conditions. Assuming that 
pairing correlations are not created during neutrino produc¬ 
tion, the amplitude is of the order of the equilibrium value, 
Eq. (88). We then find again that the mean number density 
created by pairing correlations is ~ k 2 . Therefore, the inclu¬ 
sion of the particle-antiparticle coherence leads to a negligibly 
small Sp ~ 5p ~ k 2 ~ ICG 22 . 

C. Including neutrino-neutrino interactions 

In the previous subsection we have estimated the “natural” 
size of k assuming that the neutrino-neutrino interactions are 
negligible. However, in a supernova the neutrino density is 
very large and the neutrino background may play an impor¬ 
tant role. This complicates the estimate of k because H' /!/ in 
Eq. (88) itself depends on k once we include neutrino-neutrino 
interactions, 

hT 5 = - V(3 - 2v / 2G f e* / [qf + en + eV], (91) 

J q 

see Eqs. (48) and (49). A further complication arises from the 
fact that the stationary value for k of one momentum mode p 


where we have integrated up to the cutoff |q| = Mw- Let us 
introduce the notation 

Hq P = -VP - 2v / 2G f e* f qi. (93) 

J q 

Then using Eq. (92) we can write Eq. (91) as 

H v0 « Hf ^1 - ^ . (94) 

In other words the k terms in Eq. (91) effectively lead to a 
renormalization of the total potential produced by the mat¬ 
ter and neutrino backgrounds. Numerically, the correction is 
small, \f2G-pMyv / (37t 2 ) «3x 10 -3 , and can be neglected. 

In a supernova the neutrino density is comparable to that 
of matter. Whereas each individual neutrino is relativistic, the 
bulk velocity of the neutrino background is also comparable 
to the matter velocity. Thus, the neutrino density contribution 
to Eq. (93) is not expected to be larger than the matter con¬ 
tribution. Furthermore, because the direction of the neutrino 
background flux is more likely to be parallel to the momenta 
of individual neutrinos, whereas the build up of the particle- 
antiparticle coherence requires a current component orthogo¬ 
nal to the neutrino momentum, there is an additional suppres¬ 
sion as compared to the matter effect. All in all, the estimates 
of k presented above remain essentially unaltered by the in¬ 
clusion of the neutrino-neutrino interactions. 


D. Resonance condition 

As follows from Eq. (88), particle-antiparticle coherence 
can be resonantly enhanced if H' /,y = W' ,v . In general for a 
relativistic matter flow that also includes the neutrino flux, the 
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Hamiltonian matrix reads, see Eq. (48), 


H = E 



(V° - K, —Vj_ A 
^ V° + V„) ’ 


(95) 


where, as before, V„ = pV and V± = eV are components of 
the matter flux parallel and orthogonal to the neutrino momen¬ 
tum. The resonance condition then translates into E = V ,,. 
Even assuming a relativistic matter flow, for typical supernova 
parameters V„/E ~ 10~ 9 . In other words, the resonance con¬ 
dition cannot be fulfilled in a supernova and there is no reason 
to expect k to be larger than the estimate presented above. 

Note also that for V ~ E not only does the Fermi approxi¬ 
mation break down, but also the perturbative description is no 
longer applicable. In other words it is in principle not possible 
to hit the resonance without rendering the developed formal¬ 
ism meaningless. 


E. Initial conditions 

All physical processes in which neutrinos are created have 
time scales much larger than the time scale of k oscillation. 
Hence, even during the production process neutrinos would 
adiabatically adapt to the propagation basis with respect to 
pair correlations. On the other hand, the time scales of flavor 
and helicity oscillation are much larger than those associated 
with production and detection. This separation of time scales 
is crucial for the idea that neutrinos are produced in an eigen¬ 
state of interaction, i.e., in a coherent superposition of propa¬ 
gation eigenstates. For the same physical reason, as neutrinos 
stream away from the supernova, they have enough time to 
adiabatically adapt to the external background. Thus, k does 
not oscillate but instead closely tracks its equilibrium value. 
This makes dynamical equations for k essentially superflu¬ 
ous. As the neutrinos leave the supernova, the mean pair cor¬ 
relations approach zero adiabatically and decouple from the 
evolution of p and p. 


F. Lorentz covariance 


to each other. In the frame of the first observer, the neutrino 
has momentum p along the z axis and the matter is at rest, i.e. 

V„ = V± = 0, 


H 





(96) 


Thus no particle-antiparticle coherence builds up. In the frame 
of the second observer which moves with velocity /3 along the 
x axis the Hamiltonian is no longer diagonal. 





+ V 



(97) 


and we expect helicity coherence to build up. In other words, 
physical results seem to depend on the frame. 

As we have learned from the analysis of an analogous prob¬ 
lem for helicity coherence, the kinetic equations are covariant 
only if the initial conditions also transform under the boost. 
Pair correlations “couple” neutrinos of opposite momenta. 
The notion of opposite momenta is not Lorentz invariant and 
is violated by, e.g., a boost orthogonal to the neutrino momen¬ 
tum. This alone implies that the initial conditions, which in¬ 
clude specifying k for all momentum modes, are not Lorentz 
invariant. At the same time the very fact that the definition of 
k involves two momentum modes makes it rather difficult to 
derive the corresponding Lorentz transformation rules and we 
will not attempt the derivation here. 

We have argued in the previous subsection that neutrinos 
are produced and propagate in an eigenstate with respect to 
particle-antiparticle coherence. In Sec. VI we have observed 
that if in every frame we diagonalize the Hamiltonian then the 
(transformed) eigenstate of the Hamiltonian remains invariant 
under the boosts. Here we assume that the same holds true 
also for particle-antiparticle coherence. As a particularly in¬ 
teresting example let us assume that in the frame of the first 
observer the system is in the vacuum state of the interacting 
Hamiltonian, i.e. p = p = k = 0, 


R = 



(98) 


In the early Universe, the rest frame of the plasma is the 
only natural reference frame and the question of Lorentz trans¬ 
formation properties of the pair correlators does not arise [52]. 
In a supernova environment the situation is more compli¬ 
cated. In particular, the comoving frame of the matter can 
in some cases be more convenient than the rest frame of a 
distant observer. Similarly to helicity coherence, the particle- 
antiparticle coherence builds up only if the off-diagonal com¬ 
ponents of the Hamiltonian matrix are not zero. However, 
because the off-diagonals are proportional to the component 
of the matter flow orthogonal to the neutrino momentum, their 
value depends on the frame. In particular, one can find a frame 
where the off-diagonals vanish and no particle-antiparticle co¬ 
herence builds up. 

Let us consider the same example as in Sec. VI. We have 
two identical observers moving with velocity /? with respect 


The Hamiltonian matrix Eq. (96) and the extended density 
matrix Eq. (98) commute and therefore the latter is constant 
in time. According to our assumption, after diagonalizing 
Eq. (97) by a Bogolyubov transformation, the transformed R 
takes the form Eq. (98). Transforming back to the initial basis 
we obtain. 


1 /" 1 Cfl s$ \ 

2 ( sq 1 + c $ J ’ 


(99) 


where •& is the angle of the Bogolyubov transformation that 
diagonalizes Eq. (97), tan 2$ = (J3V/'yE)/[l — /3 2 (V/E)\. 
By construction the Hamiltonian matrix Eq. (97) and the ex¬ 
tended density matrix Eq. (99) commute and the latter is time 
independent as well. 

A perplexing feature of Eq. (99) is that it seems to describe 
a state with a nonzero number of particles and antiparticles. 
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Whereas the first observer would see neither neutrinos nor an¬ 
tineutrinos, the second observer that moves with respect to the 
first one with a constant velocity /3 seems to observe a nonzero 
density of neutrinos and antineutrinos. Put in other words, the 
empty space perceived by the first observer appears to be filled 
with neutrino-antineutrino pairs in the frame of the second ob¬ 
server. However, it is not entirely clear if the (anti)particle 
densities in Eq. (99) describe electroweak interaction eigen¬ 
states and thus would actually manifest themselves via, e.g., 
particle production or momentum transfer to nuclei in scatter¬ 
ing processes. 

G. Interpretation of the Bogolyubov transformation 

To better understand the meaning of the Bogolyubov trans¬ 
formation, we solve the equation of motion for a massless 
neutrino field coupled to a constant classical current V 1 ' , and 
demonstrate that this solution reproduces the results obtained 
using the Bogolyubov transformation. 

In the Fermi limit C = i'\o f, " aa (id ll — V^)v a . Varying the 
Lagrangian with respect to the neutrino field, we obtain the 
equation of motion, a^’ aa (id^ — V t ,) v a = 0. Its solution can 
be written in a form similar to Eq. (45), 

v{t, P) =a(i,p)X-(Pv) + ^ t (* ) -p)x+(Pv), (100) 

where a(t , p) = a(t , p)e~ lul+t and (t, —p) = (—p)e*“ _t 

satisfy the usual anticommutation relations, p v is the unit 
vector in the direction of p — V, and the energy spectrum 
is given by u>± = |p — V| ± V°. 

Using the orthogonal vectors p and e we can write u>± 
in the form uj± = sj(E — pV) 2 + |eV| 2 ± V°, which repro¬ 
duces the eigenvalues of the Hamiltonian matrix Eq. (48). The 
spinor contractions, see Eq. (12), now include % (Pv)- B y 
construction, V 1 ' 1 ' and T"" vanish once we use the solution of 
the equations of motion. The diagonal elements can be ex¬ 
panded in terms of x T (p). For example for V l/V we obtain 

xl (Pv)^ M X- (Pv) = cm M (p) + Re[c 2 e /i (p)]. (101) 

Multiplied by V^, Eq. (101) reproduces the interaction part of 
the H vv element of the diagonalized Hamiltonian matrix. The 
decomposition coefficients 

E- pV , e*V 

cr = - ttt and c 2 = — -, (102) 

|p-v| |P-V| 

are related to the angle of the Bogolyubov transformation by 
Ci = cos 2$ and |c 2 | = sin 2$ respectively. In other words, 
diagonalizing the Hamiltonian matrix by a Bogolyubov trans¬ 
formation in every frame is equivalent to using the equation of 
motion. This equivalence suggests interpreting physical par¬ 
ticle densities as propagation eigenstates of the full Hamilto¬ 
nian in line with the discussion in Sec. VIIE. 

VIII. SUMMARY AND CONCLUSIONS 

Neutrino flavor conversion is important in supernovae, 
yet a full understanding remains elusive, largely because of 


neutrino-neutrino refraction and concomitant self-induced fla¬ 
vor conversion, an effect caused by run-away modes of the in¬ 
teracting neutrino gas. The difficulties in developing a robust 
phenomenological understanding of even this relatively sim¬ 
ple case explains the reluctance to add further complications. 
Yet other effects could be important as well, caused by inho¬ 
mogeneities and anisotropies of the medium and by magnetic 
fields, especially if one broadens the view to include, for ex¬ 
ample, magnetars or neutron-star mergers. It is often thought 
that helicity conversion effects will be small, at least if neu¬ 
trino dipole moments have no additional contributions beyond 
those provided by their masses, yet one should remain open to 
such possibilities. Finally, beyond flavor and helicity correla¬ 
tions, it has been stressed recently that pair correlations could 
also become important. 

Motivated by these concerns, we have studied extended ki¬ 
netic equations that describe flavor, helicity, and pair corre¬ 
lations, limiting ourselves to the mean-field level, i.e., con¬ 
sidering only propagation effects for freely streaming neutri¬ 
nos. Based on the “forward Hamiltonian” of neutrinos inter¬ 
acting with a background medium, we have derived the vari¬ 
ous terms and have given explicit results up to lowest order in 
the neutrino mass, similar to previous studies in the literature. 
For Dirac neutrinos, we confirmed previous results and have 
extended them to include magnetic-field effects. For Majo- 
rana neutrinos, we found a small correction to the mean-field 
Hamiltonian which arises from lepton-number-violating con¬ 
tractions that appear only in the Majorana case. To analyze 
the behavior of these additional terms in the limit of vanishing 
neutrino masses, we have also studied extended kinetic equa¬ 
tions for Weyl neutrinos. 

The density matrix formalism allows one to treat helicity 
oscillation induced by matter currents and by magnetic fields 
on equal footing for both Dirac and Majorana neutrinos. We 
have derived the mean-field Hamiltonian induced by electro¬ 
magnetic fields and compared it to that induced by matter cur¬ 
rents. Somewhat surprisingly, for typical supernova parame¬ 
ters, matter currents dominate over magnetic fields. In princi¬ 
ple, resonant enhancements can be achieved, for example by 
relativistic flows of matter and background neutrinos. 

Flavor and helicity oscillations can be complicated in detail, 
but they are conceptually straightforward. Their importance 
arises because charged-current interactions produce neutrinos 
in flavor eigenstates, and all interactions produce them in al¬ 
most perfect helicity states. This nonequilibrium distribution 
which is produced, for example, in the neutrino-sphere region 
of a supernova, subsequently evolves coherently and leads to 
the various flavor and helicity oscillation phenomena. 

Concerning pair correlations, the mean-field equations pro¬ 
duce similar oscillation equations. In the simplest case of 
massless neutrinos, the pair correlations are between neutrinos 
and antineutrinos of opposite momenta and the oscillations are 
between the empty state and the one filled with a neutrino and 
antineutrino. However, one probably cannot separate produc¬ 
tion from subsequent propagation. The oscillation frequency 
is here twice the neutrino energy, so in contrast to flavor and 
helicity oscillations, there is no separation of scales between 
the energy of the state and the oscillation frequency. Proba- 
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bly, as far as pair correlations are concerned, one should pic¬ 
ture neutrinos as being produced in eigenstates of propaga¬ 
tion in the medium and not as eigenstates of the interaction 
Hamiltonian. Flavor and helicity oscillations become impor¬ 
tant only because one produces a coherent superposition of 
different propagation eigenstates. As this crucial characteris¬ 
tic appears to be missing for pair correlations, we are tempted 
to suspect that pair correlations remain a small correction to 
neutrino dispersion. 

In the simplest case, helicity and pair correlations build up 
only in anisotropic media because angular-momentum con¬ 
servation forbids mixing of states with different spin. If the 
anisotropy is a convective matter current, then there is a seem¬ 
ing paradox. In the frame with the current we expect corre¬ 
lations to build up. On the other hand, we may study these 
effects in the rest frame of the medium where no correlations 
build up due to isotropy of the background. As far as helicity 
correlations are concerned, this paradox is resolved by noting 
that the handedness of massive neutrinos is not Lorentz invari¬ 
ant. Transforming both the mean-field background and the 
neutrino states to a different frame, e.g., the rest frame of the 
medium, leads to consistent physical results. For pair correla¬ 
tions, physical results must also be the same in all frames, yet 
it is less obvious how to show this point explicitly because the 
correlated modes of opposite momentum are different ones 
in every frame. Note, however, that in the supernova con¬ 
text, there is not necessarily a natural coordinate system for 
the study of neutrino propagation. Explicitly including pro¬ 
duction and detection processes, i.e., the collision terms in the 
kinetic equation, may shed more light on this question. 

The ultimate ambition of fully understanding neutrino 
propagation in dense environments and strong magnetic fields 
requires a more complete development of its theoretical un¬ 
derpinnings. Our paper is meant as a contribution toward this 
overall goal. 
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Appendix A: Chiral spinors 

Following the conventions of Ref. [62], which differ from 
the ones used in Ref. [44] by the overall sign of 7 0 , the Dirac 
matrices in the Weyl representation, which is used in this 
work, are 


= 


0 

0 


(Al) 


where a M = (l,<x) and <f M = (1 ,—cr). Here <7 is a three- 
vector Pauli matrix and 0 and 1 are 2x2 zero and unity matri¬ 
ces respectively. The chiral projectors are 

p ‘=(JS!)• p *=(S!!)- < A2 > 

The charge-conjugation matrix is 

c = -hV = ( + 0 £ ° £ ) , (A3) 

where 

'=(+1 “o 1 ) ' < A4 > 

Notice that in two-component form, one usually writes e a p 
for this antisymmetric 2x2 matrix and e a/3 for —e appearing 
in the lower right position of C. 

In the Weyl representation and with these conventions, the 
Dirac bispinors are 


v>i(p,s) 


K,sXs(p) \ 
m;,sXs(p)) ' 


Vi{p,s) = s 


f~K-sX-s( p)\ 

V K,sX-s(p) ) ’ 


(A5a) 

(A5b) 


where p is the unit vector in the direction of p, p = |p|, s = ± 
is a helicity index, and 


a r; 


p,s 


\f 


Ej - sp 
2 Ei 


6s- 



(A6) 


where Ei = ( p 2 + to?) 1 / 2 is the energy of a neutrino with 
mass rrii. 

We may describe the modes of the neutrino field in spher¬ 
ical coordinates where the momentum components are p = 
(sin 9 cos <f>, sin 9 sin <j>, cos 9). In this case, the standard two- 
component helicity spinors are explicitly 


** <f,) = ( e “ S ta|) • (A7a) 

x-(p)=(-C 4 " ! )- 


They satisfy the orthogonality condition xJ(p)X/j(p) = 6 s h- 
The matrix elements of a 1 ' are then found by direct evalua¬ 
tion to be 


xl {p)^X- (p) = n* = (1, p), (A8a) 

X+ (P)^X+ (P) = = (1, -p), (A8b) 

Xl (P)^X- (P) = -eV = -e-**(0, e), (A8c) 

xl (P)^X + (P) = -e*V* = —e^(0, e*), (A8d) 
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where e 1 ' is a polarization vector orthogonal to r>J‘. The ex¬ 
plicit components in spherical coordinates are 


e = 


/ e i(f> cos 2 § - esin 2 § 

—cos 2 | + sin 2 |) 
y — sin 6 


(A9) 


Note that the vectors rij‘ and f 1 ' depend on p, but we do not 
show this dependence explicitly to simplify the notation. 


Appendix B: Neutrino-neutrino mean-field Hamiltonian 


Because Majorana and Weyl neutrinos have two degrees 
of freedom, in many cases it is more convenient to use two- 
component notation. For Majorana neutrinos, 



and Vi 



(Bl) 


where and v\ o are two-component fields. They are re¬ 
lated by Hermitian conjugation and transform under the (^, 0) 
and (0, tj ) representations of the Lorentz group respectively. 
To emphasize the different transformation properties, the con¬ 
jugated fields, by convention, always carry a dotted spinor 
index. The spinor indices a and a are raised (lowered) us¬ 
ing the spinor metric matrices and e a ^ (e a g and e A g). 
Left-handed Weyl fields satisfy the condition P R v = v. Their 
explicit form can be obtained from Eq. (Bl) by applying the 
chiral projectors. 

Rewritten in terms of the two-component fields, the 
neutrino-neutrino Hamiltonian density of Eq. (21) is 


U vv 


Gy * 


v f • 

i,a. 


j[ 




3,p\ ‘ 


(B2) 


Taking expectation values of products of two of the four neu¬ 
trino fields and bearing in mind that fermions anticommute we 
obtain for the mean-field Hamiltonian 


^' F \ ' =-ri,aa-^/3 


IJ 


X 


2 ~ 2 V i,a V 3,P ( v lg V ^) 




+ "ta 11 ] g + v i,ocVj,P «4,d) 


(B3) 


Translating back to the four-component notation we obtain 
[Fi 7 M PLVi\(vj'yfiPLVj) for the first term in Eq. (B3). Using a 
Fierz identity [62], a^' aa a^ = — and translating 

back to the four-component notation we can represent the sec¬ 
ond term in a similar form, [vi'y 11 P R Vj\{pj')pP R Vi). By rais¬ 
ing and lowering the spinor indices and reordering the fields 


the third term can be rewritten as cr^ Q c T pv\^v\ a {vj*Vj,p)- 

Using another Fierz identity [62], = 2<5 a @8^ <*, and 

translating back to four-component notation we can rewrite 
the third term in the form 2 [vjP R Cvf](vfCP^vf). Collect¬ 
ing all terms we obtain in four-component notation 

HZ = V2G F J2(l^PLViK^PLVj) 

lJ + [Vi'l tl PL^j]{Vj'y^PL^i) 

+ [viP R CvJ](vJCP L Vi) 

+ [vf CP L Vj](iv 3 P R Cvf >) . (B4) 

Using the definition of the charge-conjugate field, v c = Cv T , 
and the resulting v c = v T C we can further simplify and write 
the last two terms in a form which coincides with Eq. (40), 

?Cf = ^G F Y^{[vi'i^P L Vi}{v j 'i ll P L v j ) 

13 +[vij fJ ’PLi'j}{vj'yi 1 PLVi) 

+ [ViP R Vj] {Vj P L vf 
+ tfPLvAfrPnv?)) . (B5) 

Thus, the effective Hamiltonian obtained using the two-com¬ 
ponent notation is identical to the one obtained using the four- 
component notation, as expected. 


Appendix C: Right-chiral currents 


For completeness, we provide the contractions of the 
Lorentz structure of right-chiral currents ( 7 ^P R ), which 
might arise in, e.g., beyond the Standard Model theories with 
right-handed currents. The contractions are 


(7^-Pr) 


VV 

ij,sh 


(7m-Pr) 


vv 

ij,sh 


( 7 ^Pr) 


w 

ij,sh 


( 7 » p r) 


w 

ij,sh 



2 V V 







0 +i<t>nH 

2p 


Tip 

n +id> rrij 
2p 



(Cla) 

(Clb) 

(Clc) 

(Cld) 
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